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ABSTRACT 

microRNAs (miRNAs) are small non-coding RNAs involved in post-transcriptional gene regulation. They 
were first identified in the free-living nematode Caenorhabditis elegans, where the miRNAs lin-4 and let-7 
were shown to be essential for regulating correct developmental progression. The sequence of let-7 was 
subsequently found to be conserved in higher organisms and changes in expression of let-7, as well as 
other miRNAs, are associated with certain cancers, indicating important regulatory roles. Some miRNAs 
have been shown to have essential functions, but the roles of many are currently unknown. With the 
increasing availability of genome sequence data, miRNAs have now been identified from a number of 
parasitic helminths, by deep sequencing of small RNA libraries and bioinformatic approaches. While 
some miRNAs are widely conserved in a range of organisms, others are helminth-specific and many 
are novel to each species. Here we review the potential roles of miRNAs in regulating helminth develop- 
ment, in interacting with the host environment and in development of drug resistance. Use of fluores- 
cently-labeled small RNAs demonstrates uptake by parasites, at least in vitro. Therefore delivery of 
miRNA inhibitors or mimics has potential to alter miRNA activity, providing a useful tool for probing 
the roles of miRNAs and suggesting novel routes to therapeutics for parasite control. 

© 2014 Published by Elsevier Ltd. on behalf of Australian Society for Parasitology Inc. This is an open 
access article under the CC BY-NC-ND license (http://creativecommons.Org/licenses/by-nc-nd/3.0/). 
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1. Introduction 

Gene expression in specific cells and tissues or at defined devel- 
opmental stages can be regulated by a number of mechanisms. 
Transcription factor (TF) binding to specific promoter sites, usually 
within the 5' upstream sequence, is currently the best understood 
mechanism of regulating gene expression. Presence of enhancer or 
suppressor DNA motifs can alter TF activity and gene expression 
under certain conditions. The accessibility of promoter regions to 
TF binding is also regulated by epigenetic modification of DNA or 
histones, involving methylation of DNA bases and/or methylation 
or acetylation of histones, resulting in changes in gene activity. Epi- 
genetic control of gene expression is increasingly studied in cancer 
development and stem cell differentiation (Tarayrah and Chen, 
2013), and an important role for DNA methylation in egg produc- 
tion and maturation of the parasitic fluke Schistosoma mansoni 
has also been demonstrated (Geyer et al., 2011). This review 
focuses on another level of gene control: post-transcriptional 
regulation mediated by microRNAs (miRNAs). miRNAs were first 
identified in Caenorhabditis elegans (Lee et al., 1993) and, with 
the advent of genome sequencing, have been identified in diverse 
organisms including humans, plants, viruses and parasitic 
helminths. While roles for some miRNAs have been elucidated, 
the functions of most are unknown. In this review we examine 
zcurrent data on miRNAs in helminths and their potential roles 
in regulating development, reproduction and survival within the 
host. We discuss the potential of miRNAs as therapeutic targets 
for parasite control and whether altered expression of specific 
miRNAs may play a role in anthelmintic resistance. 



2. Mechanism of microRNA-mediated gene regulation 

miRNAs are small non-coding RNAs of 21-25 bp in length. 
C. elegans genes lin-4 and let-7 were the first miRNAs to be identi- 
fied and characterized (Lee et al., 1993; Reinhart et al., 2000) and 
both are heterochronic genes - genes which regulate the timing 
of developmental events (Ambros and Horvitz, 1984). Worms lack- 
ing lin-4 repeat developmental events that normally occur only in 
the LI stage, in subsequent stages. Loss of another C. elegans heter- 
ochronic gene, lin-14, results in the opposite phenotype (Ambros 
and Horvitz, 1984). LIN-14 protein levels were shown to peak at 
egg hatching and decrease thereafter (Ruvkun and Giusto, 1989), 
while /in-4 showed the inverse expression pattern, lin-4 was cloned 
and found not to encode a protein, but a small RNA and the se- 
quence of lin-4 was observed to be complementary to multiple 
sites in the lin-14 3' untranslated region (UTR). From this it was 
concluded that lin-4 negatively regulates lin-14 post-transcription- 
ally, through binding sites in the 3'UTR (Lee et al., 1993; Wightman 
et al., 1993), thus paving the way to examine gene regulation by 
small RNAs. Subsequently, C. elegans let-7 was identified and 
shown to be a miRNA that regulated another heterochronic pro- 
tein-coding gene, lin-41, through binding sites located in the 3'UTR 
of the /in-41 transcript (Reinhart et al., 2000). Importantly, let-7 is 
conserved in diverse species, establishing the importance of miR- 
NAs as additional regulators of gene expression not only in worms, 
but also in higher organisms. 

In animals, miRNAs bind with partial complementarity to target 
sites that can be present in any region of an mRNA, but are pre- 
dominantly found in the 3'UTR (Battel 2009). Each miRNA can have 
many target mRNAs and a single mRNA can be regulated by multi- 
ple miRNAs. miRNAs are processed from long primaiy transcripts 
(pri-miRNA) of several hundred nucleotides, by a protein complex 
containing the RNase Drosha to yield precursor-miRNAs (pre- 
miRNA) which are stem-loop structures of around 70 nucleotides 
(Kim et al., 2009) (Fig. 1). The enzyme Dicer then cleaves the 



pre-miRNA, removing the loop region to form an imperfect duplex 
of around 22 base pairs. One strand of the duplex (that with a 
weaker base paired 5' end, referred to as the mature or guide 
strand) is incorporated into the RNA-induced silencing complex 
(RISC) and directs the complex to mRNA targets, resulting in trans- 
lational repression and/or mRNA destabilization (Chekulaeva and 
Filipowicz, 2009). For some miRNAs, the alternative strand of the 
duplex (the star strand) can be functional, but is often present in 
lower amounts or may be degraded. The specificity of miRNA 
binding is determined predominantly by complementarity of 
nucleotides 2-7 (seed sequence) of the miRNA to the target mRNA 
sequence (Bartel, 2009). As discussed below, complementarity can 
be used to predict potential miRNA targets. 



3. microRNAs in helminths 

3.1. Identifying microRNAs in parasitic helminths 

The increasing availability of genome sequence data for 
parasitic helminths has paved the way to identify miRNA se- 
quences using both computational and experimental approaches. 
The conservation of some miRNA sequences across species can 
be exploited to bioinformatically identify potential miRNAs. One 
approach taken previously by our lab (Winter et al., 2012) used 
the mature miRNA sequences present in miRBase (ftp://mir- 
base.org/pub/mirbase20) (Griffiths-Jones et al., 2008) to query 
available genome assemblies using BLASTN. The genomic regions 
flanking hits were excised and these candidate miRNA precursor 
sequences tested using RNAfold (Hofacker et al., 1994) and CID- 
miRNA (Tyagi et al., 2008). A similar approach has been automated 
by others and is known as MapMi (Guerra-Assuncao and Enright, 
2010). As well as being an effective first step in miRNA identifica- 
tion, bioinformatic searching is particularly useful in identifying 
miRNAs expressed at low levels, in restricted cell types or in devel- 
opmental stages not readily amenable to small RNA library 
preparation. 

De novo identification of miRNAs is best achieved by deep 
sequencing of small RNA libraries. Sequences are mapped to the 
genome and other RNA classes including tRNAs, rRNAs and mRNAs 
removed from the analysis. Algorithms such as miRDeep (Friedlan- 
der et al., 2012) which incorporates RNAfold, score the likelihood 
of sequences representing true miRNAs based on the folding ener- 
gies of genomic regions around the reads, similarity of these re- 
gions to miRNA hairpins and evidence of processing by Drosha 
and Dicer. In addition to miRNAs, small RNA sequencing also iden- 
tifies endogenous small inhibitory (si) RNAs and PlWI-interacting 
RNAs (piRNAs, known as 21 U RNAs in C. elegans). Each of these 
classes can be identified based on their different characteristic 
folding patterns, homology to coding sequence or for piRNAs, con- 
served 5' motif sequence (Ruby et al., 2006; Friedlander et al., 
2009). 

Small RNA sequencing projects have been carried out for a 
number of helminth species (summarized in Table 1) and small 
RNA data are likely to increase with additional genome sequencing 
projects underway. In the absence of assembled genome sequence 
for a particular species of interest, alternative approaches have 
been applied to identify miRNAs. These include mapping small 
RNA sequence reads to genomes of closely related species, as we 
demonstrated for Brugia pahangi (Winter et al., 2012), where the 
genome data from its sister species Brugia malayi were used, and 
for Echinococcus granulosus (Cucher et al., 201 1 ) where the Echino- 
coccus multilocularis genome was utilized. Identification of miRNA/ 
miRNA* duplex sequences with low RNAcofold folding energy and 
a l-2nt 3'-overhang pattern, consistent with Dicer processing, has 
also been used to identify novel miRNAs (Huang et al., 2010). Data 
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Fig. 1. microRNA biogenesis and target binding. (A) The primary miRNA transcript (pri-miRNA) is processed in the nucleus by Drosha ribonuclease. The resulting pre-miRNA 
(approximately 70 nucleotides in length) is exported to the cytoplasm where it is cleaved by Dicer to form the miRNA duplex. (B) The mature strand of the miRNA duplex is 
incorporated into the RNA-induced silencing complex (RISC) and directs binding to mRNA target sequences, often in the 3'UTR. Binding is specified by complementarity 
between the target mRNA and miRNA seed sequence (nucleotides 2-7). An example of sequence complementarity is shown in (B). 



available so far provide a platform for detailed analysis of small 
RNA expression, conservation, evolution and function. For exam- 
ple, members of the mir-36 family have been identified in Caeno- 
rhabditis species (de Wit et al„ 2009), Ascaris (Wang et al., 2011) 
and Brugia (Poole et al, 2010; Winter et al., 2012) as well as in 
platyhelminths including Schistosoma mansoni (Marco et al., 
2013) and the free-living planarian Schmidtea mediterranea 
(Palakodeti et al., 2006; Friedlander et al., 2009). mir-36 family 
members appear to be helminth-specific and are expressed only 
in adult and embryonic stages. Deletion of the mir-36 cluster in 
C. elegans results in embryonic and larval lethality, indicating an 
essential role in early development (Alvarez-Saavedra and Horvitz, 
2010). In Brugia, one highly conserved mir-36 family member was 
identified by deep sequencing and is enriched in adult female 
worms (Winter et al., 2012). Targeting this miRNA may therefore 
have potential in limiting microfilariae production. 

Comparative analysis has highlighted that, while some 
helminth miRNAs are widely conserved across diverse organisms 
and others are specific to helminths, many are novel. It is likely 
that, with additional genome data, homologues of some of these 
will be found in closely related helminths. C. elegans, C. remanei 
and C. briggsae share approximately 90% of miRNAs (de Wit 



et al., 2009). However, this number drops to only 27-30% when 
miRNAs from the necromenic nematode Pristionchus pacificus and 
gastrointestinal nematodes Haemonchus contortus and Ascaris suum 
are compared to Caenorhabditis miRNAs (de Wit et al., 2009; Wang 
et al., 201 1 ; Winter et al., 2012). The diversity of miRNA sequences 
may indicate rapid evolution of this family of small RNAs and could 
reflect roles for novel miRNAs in the more complex lifecycles of 
parasitic species and adaptation to different environments. 
Comparative analysis of miRNAs identified in S. mansoni showed 
a similar pattern: of 112 miRNAs identified, 12% were conserved 
in other flatworms and 64% were unique to S. mansoni (Marco 
et al., 2013). 

Analysis of miRNAs in P. pacificus, H. contortus, Brugia and schis- 
tosome species suggest that novel miRNAs can evolve by several 
mechanisms including gene duplication and arm switching. The 
latter results in the mature miRNA being generated from alterna- 
tive arms of the precursor (5' and 3' mature miRNAs) and examples 
have been found for changes in arm usage in different lifecycle 
stages, suggesting functional evolution (Griffiths-Jones et al., 
2011). miRNAs can be clustered in the genome and interestingly, 
clustering of some microRNAs is conserved across species, suggest- 
ing functional linkage. For example, clustering of mir-2 family 
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Table 1 

Summary of parasitic helminth species for which microRNA sequence data is 
available. 



Parasite species 



Developmental stage 



References 



Ascaris suum 

Brugia malayi 

Brugia pahangi 

Haemonchus contortus 

Trichinella spiralis 

Angiostrongylus 
cantonensis 



Bursaphelenchus 

xylophilus 
Schistosoma japonicum 



S. mansoni 



Echinococcus granulosis 

Taenia saginata 
Clonorchis sinensis 



Germline, zygote, embryo, LI , 
12 

Male and female adults 
Mixed stage (males, females, 
mf) 

Mixed sex adults, L3 

Mixed sex adults, L3 

Muscle stage larvae 

Male and female adults 

Young adults 

Mixed stages of different 
isolates 

Mixed sex adults, 
schistosomula 
Cercariae, schistsomula, 
Male and female adults 
Eggs 

Mixed sex adults 

Male and female adults 

Protoscolex 

Adult worms 
Adult worms 



Wang et al. 
(2011) 

Xu et al. (2013) 
Poole et al.(2010) 

Winter et al. 
(2012) 
Winter et al. 
(2012) 
Chen et al. 
(2011a) 
Chen et al. 
(2011b) 
Chang et al. 
(2013) 
Huang et al. 
(2010) 
Huang et al. 
(2009) 

Cai et al. (2011) 
Cai et al. (2013) 

Simoes et al. 
(2011) 
Marco et al. 
(2013) 
Cucher et al. 
(2011) 

Ai et al. (2012) 
Xu et al. (2010) 



members with mir-71 is conserved in Schistosoma species (Huang 
et al., 2009; de Souza Gomes et al., 201 1 ), S. mediterranea (Palakod- 
eti et al., 2006) and in H. contortus (Winter et al., 2012). These 
microRNAs have roles in suppressing apoptosis and in stress re- 
sponses and their conserved linkage suggests selection for their 
co-ordinated expression. 

Detailed analysis of small RNAs during early development of 
Ascaris identified some interesting features (Wang et al, 2011). 
De novo synthesis of miRNAs was detected in Ascaris embryos 
immediately following fertilization and during zygote maturation, 
a developmental stage considered to be transcriptionally quies- 
cent. These miRNAs or their precursors could not be detected in 
oocytes or spermatids, suggesting that they are not of adult origin. 
This is the earliest embryonic gene expression reported in any 
system. Surprisingly, another class of small RNAs, piRNAs, was 
not found in Ascaris nor in Brugia (Wang et al., 2011; Winter 
et al., 2012), both clade 111 nematodes, but was identified in clade 
V nematodes for which small RNA sequence is available (Caeno- 
rhabditis, Haemonchus, Pristionchus). These RNAs are important in 
spermatogenesis and silencing of mobile elements in the germline 
and it is suggested that other small RNAs, particularly 22G- and 
26G-RNAs, may carry out these key roles in Ascaris (Wang et al., 
2011) 

Study of miRNAs and piRNAs in the free-living planarian S. med- 
iterranea showed that these are upregulated in neoblasts, pluripo- 
tent stem cells which proliferate and act as progenitor cells for 
tissue regeneration (Friedlander et al., 2009). Importantly, knock- 
down of the Argonaute protein Ago-2, required for miRNA and 
siRNA biogenesis, led to a reduction in neoblast population and de- 
fects in regeneration (Rouhana et al., 2010). Surprisingly, among 
the miRNAs upregulated is let-7, which in other systems promotes 
cell differentiation rather than proliferation. In planarians, let-7 
may function differently or be involved in exiting the stem cell 
state. Recent studies identified neoblast-like cells in adult stage S. 
mansoni and it is speculated that these may contribute to the 



longevity of this parasite as well as to growth modulation of 
male-female pairs and tissue regeneration following drug 
treatment (Collins et al., 2013). Understanding the mechanisms 
regulating neoblast maintenance and differentiation and the 
potential involvement of miRNAs and piRNAs may identify new 
therapeutic targets to reduce parasite development and survival. 

3.2. Roles of microRNAs in C. elegans 

Currently the functions of only a few helminth miRNAs are 
known. C. elegans miRNAs lin-4 and let-7 were identified in genetic 
screens which demonstrated their essential roles in regulating 
hypodermal cell fates during early larval development (Lee et al., 
1993) and the transition from L4 to adult stage (Reinhart et al., 
2000). A previously unknown miRNA, lys-6, was identified follow- 
ing a genetic screen for genes regulating the asymmetrical pattern- 
ing of C. elegans ASE chemoreceptors genes (Johnston and Hobert, 
2003). However, even in C elegans, determining the functions of 
most miRNAs has proved challenging in the absence of obvious 
phenotypes following gene knockout. This is thought to be due 
to redundancy among miRNAs; approximately 60% of miRNAs 
are members of a family of 2-8 genes (Ruby et al 2006). Family 
members share the same seed sequence, although it is not known 
whether they regulate the same target genes. In addition, some 
miRNAs may function only in certain physiological or stress condi- 
tions, or may act to maintain homeostasis. It has been shown that 
levels of miRNAs from multiple C. elegans families are increased on 
infection with pathogenic bacteria (Kudlow et al., 2012). Surpris- 
ingly, reducing miRISC activity or deleting specific miRNAs en- 
hanced worm survival following infection. This suggested that 
miRNA activity, particularly in the intestine, functions to dampen 
pathogen responses that may be detrimental to worm physiology. 
This is consistent with the reduced fecundity and sterility of many 
C. elegans pathogen-resistant mutants (Miyata et al., 2008). 

Deletion mutants of 87 C. elegans miRNAs have been generated 
and clear phenotypes reported for only a small group (Miska et al., 
2007): deletion of the mir-35-41 gene cluster (mir-36 family) re- 
sulted in temperature sensitive embryonic and larval lethality, 
deletion of mir-240, mir-786 cluster increased defecation cycle 
lengths, while single deletion of mir-48 or mir-241 (both members 
of the let-7 family) affected developmental timing. A subsequent 
study examined the effects of deleting multiple members of miRNA 
families (Alvarez-Saavedra and Horvitz, 2010). This identified 
essential roles for mir-42-44 acting with other members of the 
mir-36 family, with deletion of all resulting in a penetrant embry- 
onic and larval lethal phenotype. Deletion of the mir-51 family 
(mir-51-56) resulted in defects in late embryonic development, 
particularly a lack of pharynx attachment, while essential roles in 
locomotion, growth, egg-laying and dauer formation were identi- 
fied for the mir-58 family (mir-58, 80,81,82). Members of the mir- 
58 family are orthologs of Drosophila bantam, a miRNA controlling 
cell proliferation and apoptosis. Re-introduction of individual C. 
elegans miRNAs from transgenes was sufficient to rescue the phe- 
notypes observed following loss of all family members, supporting 
the idea of functional redundancy. 

It was also suggested that miRNAs may function with unrelated 
miRNAs or with non-miRNA genes to generate robust regulatory 
networks. Brenner et al. (2010) examined the effect of loss of spe- 
cific miRNAs in a genetic background with reduced levels of all 
miRNAs, using a C. elegans mutant of argonuate-like 1 (alg-t), a 
miRNA binding component of the silencing complex. While alg-1 
mutant worms are viable, they display developmental timing de- 
fects, molting defects and early adult lethality, due mainly to 
reduction of let- 7 family activity. In this sensitized background, en- 
hanced or synthetic phenotypes were observed for 19 of 25 miRNA 
mutants tested. Brenner et al. also examined combined loss of spe- 
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cific miRNAs and hub genes, which encode chromatin-regulatory 
proteins that interact with a range of developmental signaling 
pathways. Combined knockdown resulted in sterile phenotypes 
for 4/11 miRNAs, indicating the importance of these miRNAs in 
regulating germline development. 

Data on C. elegans miRNA expression, regulation and function 
are increasing and the approaches used can be applied to parasitic 
helminths. For some parasite miRNAs, function may be speculated 
based on homology to C. elegans miRNAs, from temporal or spatial 
expression data and from target prediction programs (see below). 
However, narrowing down the list of the most likely targets and 
testing these experimentally is a major challenge. 

3.3. miRNA expression during helminth development 

Small RNA library sequence read data can be used as a first 
indicator of miRNA abundance in particular life-cycle stages (Kato 
et al., 2009). Relative expression through development or in partic- 
ular tissues can be confirmed by quantitative real-time PCR 
(qRT-PCR), miRNA microarray or Northern blot analysis. Sequence 
read data from the infective L3 stage and mixed adult stages of B. 
pahangi and H. contortus identified a number of miRNAs differen- 
tially expressed during development, while others were abun- 
dantly expressed in both stages (Winter et al., 2012). The validity 
of this approach was supported by expression data on homologues 
of C. elegans lin-4, which showed greater read numbers in the L3 
stage, and let-7, which was more abundant in adult worms, consis- 
tent with expression of their C. elegans homologues.qRT-PCR 
confirmed the temporal expression of selected miRNAs (Winter 
et al., 2012) and a detailed pattern of miRNA expression through 
development has now been obtained for 6. pahangi and H. contortus 
from recent microarray analyses (Winter et al., unpublished data). 
Arrays were designed to include all H. contortus and Brugia miRNAs 
identified from small RNA sequencing (mapping 6. pahangi small 
RNA sequences to the available B. malayi genome) and computa- 
tional analysis, as well as C. elegans miRNAs. Arrays were hybri- 
dised with small RNAs from a number of parasite stages 
including pre- and post-infective L3 larvae, L4 and adult male 
and female worms. One miRNA of particular interest is Brugia 
mir-5364, the expression of which increases significantly within 
24 hours of transfer from the mosquito to the vertebrate host 
(Winter et al., manuscript in preparation), mir-5364 is a novel 
member of the let-7 family and its rapid rise suggests a possible 
role in adaptation to the environment of the mammalian host. A 
number of interesting target genes have been predicted bioinfor- 
matically and studies are underway to experimentally test these. 

Temporal expression data on Ascaris miRNAs from RNA library 
read frequency and Northern blot analysis, showed good correla- 
tion. Of the 97 miRNAs identified, none were abundantly expressed 
throughout development (Wang et al., 201 1 ). Rather, most Ascaris 
miRNAs seem to be expressed sequentially and could be character- 
ized into four groups: enriched in germline; zygote, early and mid- 
dle embryo stages; late embryo and larval stages; and multiple 
developmental stages. 

Profiling of microRNAs through development of S. japonicum 
identified a set, including mir-t and bantam, that are significantly 
down-regulated in lung-stage schistosomula compared to cerca- 
riae, although their potential targets are currently unknown (Cai 
et al., 201 1 ). Analysis of small RNAs expressed in S. mansoni adult 
male or female worms identified 10 female enriched and 3 male 
enriched miRNAs (Marco et al., 2013). While the mir-2/mir-71 
cluster is duplicated in schistosomes, the cluster located on chro- 
mosome 5 (autosome) shows female-biased expression while that 
on the sex chromosome does not. It was suggested that gene dupli- 
cation into an autosome may overcome selection conflict on sex 
chromosomes (Gallach and Betran, 2011). 



Examination of miRNA expression patterns at the developmen- 
tal- and, where possible, tissue-specific levels may provide some 
indication as to potential biological function. As miRNAs most of- 
ten negatively regulate their target genes, an inverse correlation 
between miRNA and mRNA expression levels may act as a filtering 
step in identifying likely target genes, although it should be 
remembered that miRNAs can have a modest effect on mRNA 
levels. mRNA microarray and RNA seq data are available for a 
number of parasitic helminths and, combined with UTR annotation 
and target prediction algorithms, are valuable resources in aiding 
target identification. 



4. What are the targets of miRNA-mediated regulation? 

4.1. Identification of miRNA target transcripts 

Identification of biologically relevant targets is key to under- 
standing the functions of miRNAs. A number of computational tar- 
get prediction tools are available and are often used as a first step 
in target identification. These score the interaction of a miRNA se- 
quence with potential binding sites, usually in the 3'UTR of mRNAs. 
Scoring is based on a number of features including sequence com- 
plementarity between the miRNA seed sequence (nucleotides 2-7) 
and the target sequence, evolutionary conservation of the 
predicted miRNA binding sites, free energy of the miRNA-mRNA 
heteroduplex, stability of the surrounding mRNA and location of 
binding sites, with positioning near the start or end of the 3'UTR 
considered favorable. TargetScan, PicTar, PITA, miRanda and mir- 
WIP are some of the programs widely used, with each having 
slightly different scoring criteria (reviewed in Thomas et al. 
2010). The mirWIP algorithm was developed from verified 
miRNA-mRNA relationships in immunoprecipitation studies in C. 
elegans (Hammell et al., 2008) and therefore includes experimental 
refinement. While these programs are available for searching 
mRNAs of humans and several model organisms, including 
C. elegans, they are not yet available for parasitic species. However, 
using annotated genome data, it is possible to apply these 
algorithms to predict parasite miRNA targets. We have used this 
approach to identify potential targets of Brugia miRNAs (Winter 
et al. in preparation), exploiting the genome data available for sev- 
eral filarial nematode species (Ghedin et al., 2007; Godel et al., 
2012; Desjardins et al., 2013). 

Experimental identification of miRNA-mRNA interactions with 
components of the miRISC is increasingly used to identify in vivo 
target genes. Proteins, such as ALG (Argonaute-Like Gene) or AIN 
(ALG-1 INteracting protein) tagged with GFP can be expressed in 
C. elegans and immunoprecipitated from miRISC using anti-GFP 
antibody. mRNAs and miRNAs extracted from the precipitates are 
then sequenced or analyzed by microarray to identify enrichment 
in specific developmental stages (Zhang et al., 2009). A variation of 
this method has been developed, referred to as high throughput 
sequencing by cross-linking and immunoprecipitation (HITS-CLIP), 
in which UV irradiation is used to cross-link miRNA-mRNAs to RISC 
proteins. Unbound mRNAs are trimmed with nuclease and bound 
complexes immunoprecipitated with antibody to RISC compo- 
nents. Protected fragments are PCR amplified and sequenced to 
identify the specific miRNA binding regions in precipitated mRNAs. 
This method has been used to identify Argonaute binding sites in 
mammalian tissues (Licatalosi et al., 2008) and has been adapted 
for C. elegans (Zisoulis et al., 2010). In the absence of transfection 
methods and use of GFP pulldowns, immunoprecipitation ap- 
proaches can be adapted to identify miRNA-mRNA interactions in 
parasitic helminths using antibodies to specific miRISC proteins. 
We are currently developing this approach to identify miRNA- 
mRNA interactions in H. contortus (Marks et al., unpublished data). 
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Direct interaction between miRNA-target mRNA binding sites 
can be tested using reporter assays, measuring the level of luciferase 
or GFP expression in cells over-expressing a specific miRNA of inter- 
est. Commonly, the 3'UTR sequence is cloned downstream of the 
reporter gene and miRNA interactions ) monitored by a decrease 
in reporter activity. Candidate binding sites can be deleted or mu- 
tated and restoration of reporter levels examined. This can be carried 
out in cell culture or for nematode miRNA-mRNA interactions, in 
transgenic C. elegans. Reporter gene analysis was used to confirm 
interaction of C. elegans let-7 with the 3'UTR of Ce-lin-41, resulting 
in a decrease in lac-Z reporter activity (Reinhart et al., 2000). 

4.2. Could parasitic helminth miRNAs regulate host genes? 

miRNA target gene analyses has focused predominantly on po- 
tential targets within the same organism or tissue. However miR- 
NAs of pathogens may also have effects on host mRNA 
expression and be important in regulating host-parasite interac- 
tions. This has been demonstrated for viral miRNAs, particularly 
from human cytomegalovirus (HCMV), and may also be true of par- 
asite miRNAs. Several HCMV miRNAs are important in the estab- 
lishment and maintenance of viral infection by modulating host 
genes involved in cell cycle control (Grey et al., 2010) and in natu- 
ral killer cell function (Stern-Ginossar et al., 2007). 

The importance of miRNAs in modulating mammalian immune 
cell development and activity are increasingly understood (David- 
son-Moncada et al., 2010). For example, the miRNA signature of T 
regulatory (Treg) cells has been characterized and among the miR- 
NAs expressed are mir-21 and mir-31, which have opposing effects 
on the Treg TF FOXP3 (Rouas et al., 2009). In both B and T cells, mir- 
155 is an important regulator and loss of mir-155 can result in in- 
creased propensity for Th2 rather than Thl responses (Rodriguez 
et al., 2007). In macrophages and monocytes, mir-155 expression 
increases on stimulation with LPS and was shown to positively reg- 
ulate inflammatory responses (Tili et al., 2007). Supporting a pos- 
sible role for pathogen miRNAs in modulating host immune 
function, recent work showed that small RNAs from a common 
plant fungus, Botrytis cinerea, can silence genes involved in plant 
immunity (Weiberg et al., 2013). The fungal small RNAs bind to 
plant Argonaute proteins and mediate silencing, particularly of 
mitogen-activated protein kinase genes, to suppress immunity. 

An important finding from studies of helminth miRNAs is that 
many are novel and not conserved in C. elegans or other organisms. 
This may reflect the complex life-cycles of parasitic species, such as 
their ability to adapt to different temperatures and environments 
and longer developmental times compared to C. elegans. Interest- 
ingly, microarray analysis of H. contortus and 6. pahangi miRNAs 
indicates enrichment of novel miRNAs in the adult stage. While 
these could be involved in mating and reproduction, we speculate 
that some may act to regulate host genes to enable parasite devel- 
opment and survival within the host environment. The ability of 
parasitic helminths to suppress host immune responses is well 
documented (Hewitson et al., 2009) leading to interest in their 
use as therapeutics for autoimmune conditions such as inflamma- 
tory bowel disease, asthma and rheumatoid arthritis (Harnett et al., 
2010; Navarro et al, 2013; Weinstock and Elliott, 2013). Protein 
and glycan components of helminth excretory-secretory (ES) prod- 
ucts have been studied as potential sources of immunoregulatory 
factors. Could miRNAs present in ES or released from dying worms 
also have evolved to modulate host responses to infection? 

Surprisingly, miRNAs are extremely stable when released from 
cells and can be readily detected in serum, leading to their use as bio- 
markers for particular disease conditions including tissue injury and 
cancer (Brase et al., 2010). Recent work has also demonstrated thatS. 
mansoni-derived miRNAs in human serum have potential as bio- 
markers of infection (Hoy et al., 2014). The stability of extracellular 



miRNAs is thought to be due to their release within microvesicles or 
exosomes or by complexing with proteins, including the argonaute 
protein Ago-2, protecting them from RNase digestion (reviewed in 
Hoy and Buck (2012)). We have identified several miRNAs from 
in vitro ES material of H. contortus adult worms (Gu et al., unpub- 
lished data). One miRNA of particular interest, Hco-mir-5352, is 
upregulated in adult worms compared to larval stages and, from cur- 
rent data, homologous sequences are present only in nematodes 
infecting the gastrointestinal tract (Winter et al., 2012). The T cell 
activation marker CD69 was identified as the highest scoring poten- 
tial host target of Hco-mir-5352 and studies are currently underway 
to test for interaction between Hco-mir-5352 and CD69 3'UTR and 
any effects of this miRNA on immune cell activation and differenti- 
ation. This is potentially important for identifying novel mecha- 
nisms by which helminths may manipulate the host environment 
to aid survival, as well as identifying novel intervention strategies. 

5. miRNAs as potential therapeutic targets 

5.1. Development of miRNA inhibitors and mimics 

Levels of many miRNAs, including let-7, as well as RISC activity, 
are in general decreased in many cancer cell lines and tumours. 
miRNAs are therefore receiving increasing attention as potential 
therapeutic targets and prognostic markers for cancer treatment. 
As a result, there has been significant progress made in recent 
years in the design and effective delivery of miRNA inhibitors 
and mimics. Where available, genetic knockouts, with loss of spe- 
cific miRNAs or RISC components, or use of miRNA sponges (Ebert 
and Sharp, 2010) can also be used to determine effects of miRNA 
depletion. miRNA inhibitors, or anti-miRs, are modified anti-sense 
oligonucleotides with reverse complementary sequence to the ma- 
ture endogenous miRNA to be inhibited. Anti-miRs form a duplex 
with the endogenous miRNA, leading to it being unable to bind 
to the cellular mRNA target (reviewed in van Rooij et al., 2012). 
miRNA inhibitors can be modified to increase resistance to nucle- 
ase digestion and enhance uptake into cells. For example, inclusion 
of 2'-0-methyl (OMe) group modifications, use of locked nucleic 
acids (LNA) with bridged 2'oxygen and 4' carbon ribose moieties, 
or introduction of phosphorothioate rather than phosphodiester 
linkages between nucleotides, all enhance the stability of anti- 
miRs. Modifications are more limited with miRNA mimics, which 
are RNA duplex molecules, as the guide strand has to be identical 
to the miRNA of interest and function in the RISC. Inhibition of 
miRNA function was first demonstrated in C. elegans following 
injection of a 2'-0-methyl oligonucleotide complementary to let- 
7 into larval stages, resulting in the characteristic let-7 adult lethal 
phenotype (Hutvagner et al., 2004). Using fluorescently-labeled 
small RNAs, we have observed significant uptake following soaking 
of B. pahangi larvae and adult stages, suggesting that miRNA 
inhibition or mimickry is feasible in parasitic helminths (Fig. 2 
and Winter et al., manuscript in preparation). 

Currently, one LNA-based miRNA inhibitor (miravirsen) is in 
human clinical trials. This is an inhibitor of mammalian mir-122, 
a liver-specific miRNA which, as well as being important for liver 
function, binds to the hepatitis C virus RNA genome and is required 
for HCV replication (Jopling et al., 2005). In human patients suffer- 
ing from chronic HCV, miravirsen significantly reduced viral load, 
demonstrating the effectiveness of anti-miRs as therapeutics (Jans- 
sen et al, 2013). 

5.2. Inhibition of miRNA-mediated pathway regulation 

miRNAs are often considered to be fine tuners of gene expres- 
sion. However as demonstrated by loss of let-7 or lin-4, specific 



C. Britton et al. / International Journal for Parasitology: Drugs and Drug Resistance 4 (2014) 85-94 



91 




Fig. 2. In vitro uptake of fluorescently-labeled siRNA by Brugia pahangi. 100 mosquito-derived L3 stage larvae (A and B) or 1 adult female worm (C and D) were incubated in 
0.5 ml (for L3s) or 1 ml (for adults) of RPM1 1640/5% foetal calf serum/1% glucose at 37 °C/5% C0 2 , in the presence or absence of 0.5 uM Cy3-labeled siRNA (Ambion). After 
48 h, parasites were washed three times in 1 ml PBS/0.1% Tween and pipetted (L3s) or picked (adults) onto a 2% agarose pad on microscope slides and viewed by fluorescent 
microscopy (A and C) or under bright field (B and D) at xlO magnification. Images were collected using an Axioskop 2 Plus microscope (Zeiss), ORCA-ER digital camera 
(Hamamatsu) and Openlab (Improvision) software. No fluorescence was observed in the absence of Cy3-labeled siRNA. Scale bar = 100 pm. 



miRNAs can clearly act as developmental switches. This is ex- 
plained by the relative levels of a miRNA and target mRNA and 
the strength of their interaction, resulting in a molecular titration 
effect. For some target genes, the protein product functions only 
above a threshold level. If target abundance increases towards 
the threshold, this has the effect of titrating away the miRNAs 
available for binding, thus switching from gene repression to acti- 
vation (Mukherji et al., 201 1 ). miRNAs can act as master regulators 
of developmental transitions and responses to environmental 
changes/stresses, making them suitable control targets. Signaling 
molecules and transcription factors are particularly sensitive to 
changes in concentration and, as well as regulating expression of 
miRNAs, some TFs are also targeted by miRNAs (Leung and Sharp, 
2010). Experimental and computational prediction analysis 
showed that TF-miRNA feedback loops occur more frequently than 
expected by chance (Martinez et al., 2008). Immunoprecipitation 
studies in C. elegans also identified enrichment for receptors, ki- 
nases and ion channels among the target mRNAs in the RISC com- 
plex (Zhang et al 2009). The important roles that these dominant 
miRNA targets play in helminth development suggest that changes 
to their regulation via miRNA inhibitors or mimics is likely to be 
detrimental. 

Although miRNAs frequently exert only modest effects on target 
mRNA levels (average of 2-fold decrease), because they often reg- 
ulate many genes functioning together in the same pathway or 
protein complex, their effect can be amplified. In addition, the ef- 
fect of miRNAs on target gene level can be much greater in specific 
cells where the abundance of some miRNAs is high. With increas- 
ing data on temporal miRNA expression, as well as identification of 
many helminth-specific miRNAs, we can begin to test inhibitors 
and mimics as novel therapeutics. Will inhibiting parasite miRNAs 
upregulated on host infection, such as Bpa-mir-5364, or in adult 
worms limit parasite survival and reproduction? Not only are miR- 
NAs potentially valid targets, but understanding the overall effects 



of miRNAs by characterizing their mRNA targets and their tempo- 
ral and spatial expression, should identify co-ordinated pathways 
and networks required for developmental transitions or responses 
to stresses, which may themselves be targets for parasite control. 



5.3. Involvement of host miRNAs in parasitic infections 

Not only is understanding roles of parasite miRNAs important 
to infection control, but studying host miRNAs may also lead to no- 
vel intervention strategies. Characterizing changes in host cell 
miRNA profile following parasite infection may elucidate triggers 
and mechanisms regulating host immunity, pathology and inflam- 
mation as well as parasite-induced suppression. There is increasing 
evidence that miRNAs can regulate Toll-like receptor (TLR) signal- 
ing molecules, such as AP-1 and NFk-b, as well as TLR expression. It 
is speculated that miRNAs may provide a link between innate and 
adaptive immune pathways (reviewed in O'Neill et al., 2011) and 
are therefore relevant to understanding regulation of host immu- 
nity. A number of studies have examined changes in host miRNA 
levels following infection to identify parasite-induced changes 
and regulation of host responses. For example, infection of leuko- 
cytes with the transforming parasite Theileria results in an increase 
in host mir-155 which, via TFs AP-1 and c-Jun, acts to maintain the 
cell hyperproliferation state and parasite survival (Marsolier et al., 
2013). Inhibition of mir-155 or AP-1 reverses parasite-induced 
transformation. Changes in miRNA expression also occur in host 
macrophages following implantation of 6. malayi adult worms into 
the peritoneal cavity and this has been shown to regulate alterna- 
tive macrophage proliferation (Ruckerl et al., 2012). Recently it was 
shown that invasion and induction of innate immunity in myeloid 
cells by North American eastern equine encephalitis virus is 
blocked by a haematopoietic-cell-specific miRNA (mir-142-3p). 
This host miRNA binds to the 3' end of the viral RNA genome to 
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restrict its replication in this cell lineage, thus suppressing the im- 
mune response and aiding virus survival (Trobaugh et al., 2013). 

6. microRNAs in drug resistance 

miRNA expression or activity are often altered in cancer cells 
and some of these changes are also thought to be associated with 
development of drug resistance. Drug treatment of cancer cells 
has been shown to result in altered expression of miRNAs which 
target drug metabolizing enzymes and drug transporters leading 
to drug resistance (Yu, 2009). For example, gemcitabine treatment 
of xenograft tumours results in decreased expression of mir-328, a 
miRNA which negatively regulates the ABCG2 efflux transporter 
(Pan et al., 2009). Decreased mir-328 levels therefore lead to in- 
creased transporter expression and multidrug resistance. 

We have speculated that the link between changes in miRNA 
levels and drug resistance in cancer cells may also be a feature of 
drug resistance in parasitic nematodes (Devaney et al., 2010). 
Altered levels of drug targets, including ion channels and nAChR 
subunits have been found in isolates of H. contortus resistant to 
ivermectin (Rao et al., 2009) and in A caninum resistant to pyrantel 
(Kopp et al., 2009). Following drug selection, changes were also 
found in expression level and allele frequency of ABC transporters 
(Prichard and Roulet, 2007) and P-glycoproteins (Blackhall et al, 
2008). Whether these changes in gene expression are mediated 
by miRNAs is currently unknown. In C. elegans, mutants of mir-1, 
which negatively regulates nAChR subunits unc-29 and unc-63, 
show decreased muscle sensitivity to Ach and resistance to levam- 
isole (Simon et al., 2008). This is thought to result from changes in 
subunit composition affecting the biogenesis or function of hetero- 
multimeric receptors. Recently, we have used miRNA microarray 
analysis to examine the expression profile of C. elegans strains 
and H. contortus isolates resistant to ivermectin. Our initial data 
have identified significant and consistent changes in expression 
of three miRNAs - Cel-mir-85, Cel-mir-788 and a novel H. contortus 
sequence Hco-mir-9551 (Gillan et al., unpublished data). Work is 
underway to identify alterations in mRNA levels which may be 
associated with these miRNA changes. Obvious targets include 
drug receptors, transporters and drug metabolizing enzymes, but 
the transcriptomic approach has the potential to identify novel 
mechanisms of drug resistance. Alternatively, the observed 
changes in miRNA expression levels may be a consequence of other 
changes associated with resistance. These miRNAs may also prove 
to be valuable and sensitive markers of anthelmintic resistance. 

7. Conclusion 

microRNAs of parasitic helminths are receiving increasing 
attention as novel regulators of development, sex differentiation, 
drug resistance and host-parasite interactions. The important and 
varied roles of miRNAs identify these as valid targets of interven- 
tion strategies. With genome sequence data available for more par- 
asite species, miRNA identification and understanding of miRNA 
evolution and function will increase. Which miRNAs are conserved 
across species, which are specific to particular parasite clades or to 
modes of transmission or location within the host? Improved gen- 
ome assembly/annotation together with comparative analyses will 
also aid target gene predictions through identification of conserved 
miRNA binding sites. Parasite miRNAs may also regulate host tar- 
get genes and there is increasing evidence that pathogens exploit 
host miRNAs to aid replication and survival. This has been demon- 
strated for HCV, which relies on host mir-122 in liver cells, Theileria, 
which alters leukocyte mir-1 55 levels, and equine encephalitis 
virus where cell tropism seems to be mediated by a host miRNA. 
It will be important to examine effects of helminth miRNAs on both 



parasite and host gene expression, any parasite-induced changes in 
host miRNA levels and whether host miRNA expression could also 
influence helminth tropism. This has potential to lead to therapeu- 
tics targeting parasite pathways as well as altering host-parasite 
interactions that influence parasite survival. 
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